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SUMMARY 


The Denavit-Hartenberg parameters define the relative location of joint axis 
systems in a robot arm. A recent criticism is that one of these parameters ap- 
proaches infinity as the rotational axes of two successive joints approach a parallel 
condition. This causes an ill-conditioned transformation matrix and locates a joint 
axis system far. away from the joint itself. This paper introduces a simple modifica- 
tion of these parameters to easily overcome this criticism. This modification (which 
entails the constraint that a transverse vector between successive joint rotational 
axes be normal to one of the rotational axes, instead of both) leads to modified 
Denavit-Hartenberg parameters that more favorably locate successive joint axis 
systems . 

An example is given with respect to the rotational axes of the elbow and shoul- 
der joints in a robot arm. The regular and modified Denavit-Hartenberg parameters 
that locate the axis system of the elbow joint relative to the axis system of the 
shoulder joint are extracted in the example by an algebraic method with simulated 
measurements of three different locations of a point on a robot arm. For small mis- 
alignments of the elbow- and shoulder- joint rotational axes from a parallel condi- 
tion, the regular Denavit-Hartenberg parameters, unlike the modified parameters, were 
found to be extremely sensitive to measurement accuracy. 


INTRODUCTION 

The motion of a robot hand is the result of joint movements in the robot arm. 

To transform operator commands to the robot hand into joint movements and to pass 
sensor information along the arm, the relative location of successive joint axis sys- 
tems must be known. By far, the most popular way to describe the relative location 
and orientation of one joint axis system with respect to another is to use the well- 
known Denavit-Hartenberg parameters (ref. 1), which define the elements in transfor- 
mation matrices. However, these parameters have been criticized recently (refs. 2 
and 3) as unsuitable in the case where successive rotational axes approach a parallel 
condition. Specifically, several of the elements in the transformation matrix 
approach infinity, and a joint axis system is located far away from the joint itself. 

The purpose of the present paper is to modify the Denavit-Hartenberg parameters 
to overcome the aforementioned criticism. The algebraic method in reference 4 is 
extended to extract the modified parameters. 


SYMBOLS 



measurement vectors from world coordinate system to a point 

on 

robot arm 

A i-1 

homogeneous transformation matrix 

from coordinate system i 

to 

i - 1 

a i 

length of itj_ 




a i 

common normal vector between Z^_-| 

and Z^; intersects Zj_ 

at 

°i 



C. vector from world coordinate system to center of circular trajectory of a 

point on robot arm about line of rotation of joint i 

D. constrained transverse vector from Z. - to Z- 

i i-l 1 

E elbow of robot arm 

e^ unit vector normal to both Zj^ and Z^ that, by definition, points 

along X^ 

H hand of robot arm 

location of 0^ from in coordinate system i - 1 

i integer 

vector from world coordinate system to a point on line of rotation for 
joint i 

N neck of robot arm 

N vector normal to circular trajectory plane of point on robot arm 

0 base of robot arm 

0^,0 W origin of coordinate system i for joint i + 1 and of world coordinate 

system, respectively 

Q, Q-j t Q 2 'Q( x ' Y/ z ) points in three-dimensional space 

R. vector from 0 TT to 0. 

x w 1 

rotational transformation matrix from coordinate system i to i - 1 
r^ coordinate along Z^_-j of 0^ 

r^ vector from O.^ to coordinate along Z x-1 of °i 

S shoulder of robot arm 

n t unit vector along positive rotational axis of joint i 

normal vector between lines of rotation for joints 2 and 3 
W wrist of robot arm 

X^ axis directed along common normal between and Z^ 

X W ,Y W ,Z W world coordinate axes 

Y^ axis directed to complete right-hand axis system with X^ and Z^ 

Z^ axis of positive rotation of joint i + 1 

angle between Z^_-j and Z^, measured positive about X^ 

3^ constant bias angle, which yields joint angle 0! when summed with joint 

angle 0^ 1 
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->■ 4- “V 

joint angle corresponding to measurement vectors A, B, and C, 
respectively 

joint angle with initial value (0°) corresponding to initial position of 
robot arm 

joint angle between Xj^ and X^, measured positively about 
incremental positive change in 0^ 
real variable in vector line equation for joint i 

radius of circular trajectory of point on robot arm about line of rotation 
of joint i 

components of D. with respect to coordinate system i - 1 when 

^ = o° 1 

Mathematical notations: 

• dot or scalar product 

x cross or vector product 

T superscript to indicate transpose 

| j length of a vector or absolute value 

ANALYSIS 

The regular Denavit-Hartenberg parameters are first described and then the dif- 
ficulties encountered with their use in practical robotic applications are discussed. 
The basic set of parameters is then modified such that the criticism is no longer 
valid. Next, an algebraic approach, adapted from reference 4, is extended to extract 
the modified set of parameters. Parameters are then extracted for near-parallel 
rotational axes in a simulated robot arm. 


6 A' 9 B' 0 C 

e i 

0 

1 

A0 . 

1 

x i 

Pi 


Denavit-Hartenberg Parameters 

In figure 1, a^ (component of a^ along X^), r^ (component of r^ along 

z i-l)' a i' and referred to as the Denavit-Hartenberg parameters, completely 

characterize the relative location of successive joint axis systems. 

Joint axis systems .- Figure 1 illustrates the axis systems associated with 
joints i and i + 1. By convention, joint i is associated with the coordinate 
system i - 1. Hence, joint i corresponds to the axis system with origin at 0 i-1 , 
whereas joint i + 1 corresponds to the adjoining axis system with origin 0^. By 
definition, the axis of rotation for joint i always lies along the associated 
Z i-1 . The vector , directed toward Z^, is the normal vector between and 

Z^. The intersection point of with Z^ locates the origin 0^. The axis X^ 

originates from 0^ in the same direction as a^. In the event that Z^_<j anc ^ 

Z^ intersect (fig. 1(b)), a^ is the zero vector, and Xj_ is then directed from 
this intersection in the direction of the cross product obtained by multiplying a 

unit vector along by a vector along Zj_« The vector r^ is from the 
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origin 0 i _ 1 to the intersection of with Z i _ 1 (fig. 1(a)); for intersecting 

lines of rotation, r^ is a vector along from 0 i _ 1 to 0 i (fig. 1(b)). 

The angle is the angle between a line parallel to through the origin 0- 

and Z ir being measured positive about (fig. 1). Finally, the joint angle 0.' 1 

is the angle between X^_.| and a line parallel to through 0^_1 and is mea- 1 

sured positive about Z i _ 1 (fig. 1). For clarity, the axes and Y i _ 1 , which 

simply complete right-hand axis systems, are omitted. 

Relationship between joint angles ^ and 9? .- In general, the joint angle 9.' 
is not equal to the joint angle 0 i , which is referenced to the initial position of 
the robot arm (e.g., the position of the robot arm in fig. 2). Corresponding values 
of 0£ and 0!^ are assumed to be related by the linear equation 


0 ! 

x 


e. + e. 

X X 


( 1 ) 


where 0^ is a constant bias, reflecting an initial offset in 0!, because 0.' and 
0^ are generally measured from different starting positions. The joint angle 1 0- 
is measurable prior to establishing the joint axis systems, but 0! is not yet 
measurable. 1 

Basic coordinate transformation .- The relative joint geometry dictates the basic 
transformation equations between adjacent joints. The coordinates of an arbitrary 
point Q(x,y,z) with respect to the coordinate system i can be transformed to co- 
ordinates of Q with respect to the coordinate system i - 1 by using the relation 


with 


where 


"X" 


~x' 

y 

i 

y 


= A - , 


z 

1-1 

z 

.1. 

i-1 

_1 . 


( 2 ) 


*i-1 




■f 

h. 


(3) 


R i-i 


cos 0! 

l 


sin 0! 

l 


-cos a. sin 0.' 
1 1 


cos a- cos 0! 
i i 


sm a. 


sin a- sin 0! 
1 i 


-sin a- cos 0? 
i i 


cos a. 


(4) 
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to i 


1 


is the rotational transformation matrix that relates coordinate system i 
with respect to orientation, and where 


*>■ 




cos 0.*'' 

l 


sin 0! 

i 


(5) 


is the location of Ch from Ch * in coordinate system i - 1. As can be seen, the 

1 -L— l 

elements in the basic transformation matrix A^ are expressed in terms of the 
Denavit-Hartenberg parameters. 


Criticism of Denavit-Hartenberg Parameters 

A criticism of the Denavit-Hartenberg parameters is that | r^ | -► 00 as + 0°. 

(See fig. 1(b).) Indeed, this criticism is justified in that this behavior causes: 
(1) Sensitivity of the parameter r^ to errors in misalignment from a parallel con- 
dition (a^ = 0°); (2) Ill-conditioned transformation matrices ( r ^ terms and matrix 

products approach infinity as ■* 0°); and (3) Excessive displacement of the axis 
system from the robot arm. The intent in this paper is to eliminate these weaknesses 
by modifying the Denavit-Hartenberg parameters. 


Modified Denavit-Hartenberg Parameters 

In this paper, the regular Denavit-Hartenberg parameters are modified by only 
insisting that a transverse vector between successive joint rotational axes be normal 
to one of the axes instead of to both axes. This simple modification leads to a more 
favorable location of successive joint axis systems. 

As with the regular Denavit-Hartenberg parameters, let be the axis of rota- 

tion for joint i + 1, let X^ point in a direction that is normal to both 
and Z^, and let complete the right-hand coordinate system. The angle between 

Z^_ 1 and Z^ is a^, measured positive for rotation about . The difference 

between the regular and modified parameters is in how 0^ is located from • 

In figure 3, vu and u^ +1 are unit vectors along the lines of rotation for 
joints i and i + 1, respectively. A^so in figure 3, (determined later) is 

a vector from 0 W to a point oi£ Z^; and locate and Cb from 

0 , respectively. The vector D^ is from a point on to a point on Z^. The 

origin 0^ is removed from by 


-»* 



•> 



+ 


4 * 



( 6 ) 
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where 


D i = R i - R i _ 1 - r. Ui 


(7) 


and where r^ is the coordinate of 0^ along • 

When used in equation (3), h^ is expressed in the coordinate system i - 1. 
When 0 ± = 0°, 


/ { * 


h. 


w 


■i 

+ r, 


(8) 


where (using and e^_^ for the condition 0^ = 


0 °) 


•i = D i * e i-i 


(9) 


n i = °i * (3. X l i _ 1 ) 


( 10 ) 


and 


r. = d. • u. 

i i 


( 11 ) 


denote the projections of along X i-1 f , and 

is a unit vector along . The parameters r\^, 

place of the regular Denavit-Hartenberg parameters a^ 
rotates, the location of from 0.^ changes with 


2.^ , respectively, and e^ 
and r^ now take the 
and r^ . As joint i 
0j_ as 



cos 

0 i 

- 

q 

sin 

■* 

f 





h i = 

5^ sin 
1 1 

CD 

H* 

+ 

q 

cos 



q 

+ 

r i 



e A 

6 i 

/ 


( 12 ) 


which is used in equation (3). Equation (12) reduces to equation (5) if %. = a if 
hi = q = 0, and 0 i = 0£. 1 
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Transverse vector normal to both Z^ and Z.^ for regular Denavit-Hartenberg 

parameters When constrained to be normal to both Z i _ 1 and Z ir D i is located 
from Q-j to Q 2 in figure 3. Then, r^ becomes the location of the point Q-, 
along • The direction of is along Xj_, and its component is a^. Since 

X^ is simply rotated by 0! about Z^_.j , with components in coordinate 

system i - 1 is 






cos 0! 
1 


sin 0! 
1 


(13) 


where 0! is 0. measured from an initial zero position where X^_.| and X^ are 

parallel. Then, from equation (6), h^ is the same as equation (5) for the regular 
Denavit-Hartenberg parameters. 


The problem with constraining to be normal to both Z^_ 1 and Z^ is that 

r^ can shift 0^ far away from its associated joint on the robot arm. Conse- 
quently, two alternate constraints are imposed on in this paper. They are 

(1) Only constraining to be normal to Z^_ 1 and (2) Only constraining to 

be normal to Z^. 

*»■ 

Transverse vector normal to Z.^ .- The transverse vector is given by 

equation (7), with 


R i " C i+1 + X i+1 u i+1 


(14) 


where C i+1 (discussed subsequently) is a point on Z^. The value of A i+1 , which 
simply moves the terminal point of along Z^, is chosen to satisfy the 

constraint 


-► 




(15) 


That is. 


D i 


is normal to Z^^^ 


Equations (7), (14), and (15) show that 


A 


i+1 



■* 

+ r. u. 
1 1 


-► 


U i+1 



(16) 


Equations (9) to (12) then yield h^ . Since is normal to Z^_ 1 , = 0. The 

parameter r^ is now arbitrary and is chosen to shift 0^ to a favorable location. 
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Transverse vector normal to Z^.- The transverse vector is given by equa- 
tion (7), with given by equation (14), where is chosen to satisfy the 

constraint 


+ 

D* 


5 i + i - 0 


(17) 


That is, is normal to Z^. Equations (7), (14), and (17) yield 


X i+1 = (*i-1 + r A - C i+1 ) • 


i+1 


(18) 


Equations (9) to (12) then yield h^, and r^ is a free parameter that is chosen to 


shift 0^ to a favorable location. 


Algebraic Method To Extract Modified Denavit-Hartenberg Parameters 

Measurement data needed to extract the relative joint parameters is generated by 
individually varying the joint angles in the robot arm and measuring the location of 
a point on the robot hand or other extension. (See ref. 4.) Each joint angle is 
measured relative to a selected zero position; for example, see the position of the 
robot arm in figure 2. 

Any point on the robot arm (except a point on the line of rotation) generates a 
circular trajectory as the joint angle 0^ is varied, and the location of a point on 
the robot arm is measured relative to a world reference axis system. As subsequently 
shown, three locations of the point on the circular trajectory, along with the corre- 
sponding values of 0^, are enough to determine the unit vector u^ along the line 
of rotation for joint i and the vector from the world axis system to the cen- 
ter of the circular trajectory of the point. With and for successive 

joints, the relative joint parameters are computed. First, joint i is rotated to 
obtain u^ and C^. Then, with 0^ = 0°, joint i + 1 is rotated to obtain 
and • The process is initiated by assuming a fictitious joint with a rotational 

axis along Z w . 

Circular trajectory center .- Let A, B, and C denote three positions of a 
point on a robot arm as the point moves in a circular trajectory about the line of 
rotation for joint i (fig. 4), and let 0 A , 0 B , and 0 C be the corresponding val- 
ues of 0^. The objective is to compute the vector from the world axis system 

to the center of the circular trajectory. Let the radius of the circle be p^, then 






2 

p i 


(19) 






( 20 ) 
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(i - i t ) . (f - = 0 1 


( 21 ) 


That is, all points are located at a distance from the center of the circle. 

Equations (19) to (21) constitute three equations in four unknowns ( and the three 
components of C. ) . A fourth equation (mistakenly neglected in ref. 4) results from 

1 4 

the constraint that the tip of lies in the plane of the data points; that is, 

one of the equations 


/ 4 4 % 4 

(c^ - AJ • N = 0 

/ 4 4 

(c i - bJ • n = o 

-y 

[C i - Cj • N = 0 


( 22 ) 


(23) 


(24) 


where 


4 - . 4 - 4 - , 4 - 4 

N = (C - A) x (B - A) 


(25) 


is a normal vector to the circular trajectory plane. 

Eliminating p^ from equations (19) and (20) by using equation (21) results in 
the following equations: 


(A-C.) • (A - C,) = (S-S.) • (C-C.) 
(B-^) - (B-$.)= (C-£.) • (*-$,) 


(26) 


(27) 


Equations (22), (26), and (27) are three linear equations in three unknowns (the com- 

4 * 4 

ponents of C^) • An equivalent matrix equation for computing the components of 
is 



N ( 1 ) 


N( 2 ) 

N(3) 





i 

4- 4 

N • A 

\ 


c(i) - 

A( 1 ) 

C( 2) - A( 2 ) 

C( 3 ) - 

A( 3 ) 


C i (2) 

= 

& 

4* 4* 

• C - A 

• A) 


C ( 1 ) - 

B( 1 ) 

C( 2 ) - B( 2 ) 

C( 3 ) - 

B( 3 ) 


l C i (3) i 



4- 4 

* C - B 

• B) j 

J 

Once 

c. is 

known , 

any one of 

equations 

(19) 

to (21) 

yields 

2 

p i- 



(28) 
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Unit vector Figure 5 shows the circular trajectory of a point on the robot 

arm and two position vectors A and B, along with the incremental joint angle 


A0i = e B - 0 A 


(29) 


between these position vectors. A unit vector normal to the plane of the circular 
trajectory and passing through point (whose coordinates are the components of 
the vector C^) is 


* (A - C.) x (*- £.) 

u i = 2 

p. sin A0. 


(30) 


With 0° < A8^ < tt, u ± in equation (30) is in the same direction as the rotational 

vector of joint i. An average u^ should be used to reduce the effects of measure- 

ment errors. 

Axis X^.- The direction of X^ is defined by a unit vector e^ (fig. 3) that 

is normal to both the line of rotation for joint i and the line of rotation for 
joint i + 1 • Such a unit vector is given by 

x " i+1 )/K * " i+1 | < 3i > 

if u. x * 0, or by 

«i = Di/I^il ( 32 ) 

if u^ +1 * u^ = 0. The vector cross product in equation (31) gives a vector which is 

normal to both u^ and u^ +1 if the lines of rotation for joints i and i + 1 

are not parallel. If the rotational lines are parallel, equation (32) is used. 

Alternatively, e^ is any nonzero solution to the simultaneous linear 
homogeneous equations 


+ 




(33) 


and 


i+1 



( 34 ) 


whether or not 


+ 



and 


i+1 


are parallel. 
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Axes Z i and Y^- By convention, Z i points in the direction of positive 
rotation for joint i + 1, that is, in the direction of u^ +1 • Y^ completes a 
right-hand axis system with and Z^ . 

Equation for tan 0|.- The equation for tan 6^ is 


tan 0 ! 

i 



-»■ 

e i-1 






( 35 ) 


The numerator in equation (35) shows the cross product of a vector e^ - along X^_^ 
and a vector along X^. Forming the dot product of this result and a unit 

vector u-_ along Z i-1 produces sin 0! with the correct sign for a positive rota- 
tion about Z i _ 1 (or equivalently . The denominator is equivalent to cos 0^ . 

Hence, the fraction represents tan 0^, where 0° < 0^? < 2 tt. The joint angle 0^ in 

equation (35) corresponds to the fixed position of joint i after u^ has been 
determined and when joint i + 1 is being varied to obtain u^ + .j • The bias 0^ in 
equation (1) is just the calculated value of 0^ when 0^ = 0°. 

Tan with in direction of e^.- The appropriate equation is 


tan cl 


[ ( 5 i x ^i+i ) 



u i+i 


(36) 


The right-hand side of this equation shows the cross product of a vector along ^j __1 
(or £.) and a vector along (or u • ) , and then the dot product of this result 

and a unit vector along X£ (or ep gives sin a^. The dot product in the denomi- 
nator yields cos a^. 

Parameters , n. , C, , and r . At 6, = 0°, the components of d. in coor- 

XIX X 

dinate system i — 1 are T}., and £. in equations (9) to (11), respectively, 

and r- is a free parameter 1 chosen for favorable location of 0. . The transverse 

vector is computed in equation (7) with from equation (14), where 

depends on how is constrained. If D i is constrained normal to Z^ , equa- 
tion (16) is used for / but is constrained normal to Z i _ 1 , equa- 
tion (18) is used. In use, the constant parameters £. , h . , C./ and r appear in 

iii l 

equation (12) for h- , which only varies with the joint angle 0. (or 0* - 0., 
eq. (1)). The vector h^ is then used in equation (3). 


EXAMPLE 

For the robot arm depicted in figure 2, this example focuses on determining the 
relative location of the elbow-joint axis system with respect to the shoulder- j oint 
axis system. A transverse vector between the rotational axis 7*^ of the elbow joint 
and the rotational axis Z 1 of the shoulder joint is used in locating the elbow- 
joint axis system. 
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Transverse-Vector Constraints 


Three ways to locate the axis system result from three separate constraints on 
the transverse vector, which is constrained to be either (1) Normal to both rota- 
tional axes Zj and Z 2 , (2) Normal to the elbow— joint rotational axis Z 2 , or 
(3) Normal to the shoulder- joint rotational axis Z 1 . The first constraint leads to 
the well-known Denavit-Hartenberg parameters. 

Transverse vector normal to bo th Z 1 and Z 2> - The rotational axes of the 

shoulder and elbow joints in figure 2 are parallel. The Denavit-Hartenberg parame- 
ters are c » 2 = 0°, r 2 = SN along Z-j , a 2 = ES in the X 2 direction, which is 
normal to both Z^ and Z 2 , and 0 ^ = 0 £ + S 2 , where 6 2 = 90° is the angle be- 
tween X 1 and X 2 • These parameters locate the elbow— joint axis system as shown at 
the point E. Such is not the case, however, when the axes are not parallel. For 
example, suppose Z 2 is rotated counterclockwise by 0.1° about Y 2 . The Denavit- 
Hartenberg parameters then become a 2 = 0.1° (angle between Z y and Z 2 ), 
r 2 = -9734 in. (intersection of Z 2 with Z-, along Z 1 ) , and a 2 = 0 (since Z 2 
intersects Z 1 ) . Thus, the elbow-joint axis system is moved away (by r 2 ) from the 
robot arm by about 81 1 ft. This does not happen with the modified parameters intro- 
duced in this paper. Indeed, the elbow-joint axis system remains located at the 
desired point E. 

Transverse vector normal to Z 1 .- When the transverse vector is normal to Z^ 

(fig. 6 (a)), the point E is located by distances 5 2 along X 1 , n 2 along Y 1 , 
an< ^ r 2 ^2 along Z^ , where r 2 is now the desired location of the elbow-joint 
axis system. For this example, in which a 2 = 0.1°, these parameters are as follows 
?2 = h 2 = ES = 17 in. , ? 2 = 0 (by the constraint on the transverse vector, g 2 

is always zero and need not be determined), and r 2 = SN = 6 in. (by choice). ^ 

Transverse vector normal to Z 2 .- When the transverse vector is normal to Z 2 

(fig. 6 (b)), the parameters that locate the elbow-joint axis system are as follows: 
? 2 = °> h 2 = ES = 17 in., ? 2 = -ES/tan a 2 = -0.030, and r 2 = SN -£ 2 = 6.060 in. 
Hence, r 2 + ? 2 = SN = 6 in. 


Simulated Measurements 

In figure 7, assume that the waist- and shoulder- joint axis systems have been 
determined. Now the problem is to determine the elbow- joint axis system by using 
simulated locations of a point on the robot arm for three distinct elbow-joint 
angles. Specifically, locations of the point W are assumed to be measured with 
respect to the waist-joint axis system (which, in this example, coincides with the 
world reference axis system) for three distinct values of 63 . The simulated mea- 
surements are based on the following dimensions (in inches): NO = 26, SN = 6 , and 

ES = WE = 17. Examples of simulated measurements are as follows: 


( 0.000 
6.000 
60.000 
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These measurements give the location of point W to the nearest one-thousandth of an 
inch as 63 takes on the values of 0 A , 0 B , and 0 C . In computing the simulated 

vector positions A, B, and C (see appendix), the fact that <x 2 = 0,1° is assumed 
to be unknown. The following vectors are known: 




These vectors are expressed in the waist (or world) axis system. The vector from 0 
to N is R 1 ; c 2 is the vector from the center of the circular trajectory of point 
W as the shoulder joint is rotated; e.j is a unit vector normal to Zq and 
that defines X 1 ? and is a unit vector in the direction of rotation of joint 2 

(shoulder) in base coordinates. 
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Also, based on the simulated measurement data 


and 



u 3 



where c 2 is a vector to the center of the circular trajectory from the waist coor- 
dinate system and u^ is a unit vector along (the line of rotation of the elbow 

joint) in the waist coordinate system. 


Location of Origin of Elbow-Joint Axis System 
by Denavit-Hartenberg Parameters 

^ For the regular Denavit-Hartenberg parameters, the transverse vector, which is 
D 2 in the notation of the text, is normal to both of the lines of rotation of the 
elbow^and shoulder joints. Specifically (from ref. 4, pages 9 to 1 1 , with i = 2 
and d 2 = v*), 


D 2 = - t 2 (37) 


where 





C 3 + 


X 3*3 


(38) 


and 


-► 






X 2”2 


(39) 


with 


& - cj . [5, 


i - (1 





( 40 ) 


1 4 



(41) 


X 2 


- Cj - (S, - (u • L& 


~ & ‘ u J 


That is, A 1 and A 2 are such that d 2 • u 2 = 0 and d 2 • u 3 = 0. 

The regular Denavit-Hartenberg parameter r 2 is the component of 


r 2 ^2 ~ R 1 


(42) 


along Z 1 ; that is. 


r 2 = r 2 * u 2 


(43) 


The regular Denavit-Hartenberg parameter a 2 is the component of 


a 2 = ^3 ~ A 2 


along X 2 , that is. 


a 2 = a 2 * e, 


(44) 


(45) 


Based on the simulated measurement data, r 2 = -9789.719 in. and 
a 2 = 0.298 in. The true values are r = -9734.273 in. and a 2 = 0 in. Hence, the 
calculated values of r 2 and a 2 are in absolute error by 55.446 in. and 0.298 in., 
respectively. 

Different values of a 2 are shown in the first column of table I, and the 
Denavit-Hartenberg parameters r 2 and a 2 that locate the origin of the elbow- joint 
axis system are shown in the second and third columns, respectively. Note that 
a„ = 0 because the lines of rotation of the elbow and shoulder joints intersect. 

For measurements (components of A, B, and C) rounded off to the nearest ten- 
thousandth, one-thousandth, and one-hundredth of an inch, table I shows the absolute 
errors in the computed values of the parameters r 2 and a 2 * e lbow and 

shoulder joints approach a parallel condition (a 2 0 °), two important observations 

from table I are as follows: (1) The elbow-joint axis system is located far off the 

robot arm (large value of r 2 ), and (2) The error in the computed value of r 2 is 
excessively large, even for very accurate measurements. 
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Location of Origin of Elbow-Joint Axis System by Modified Parameters 

In this example for the modified parameters, the transverse vector from the 
elbow- joint axis of rotation to the shoulder- joint axis of rotation is normal to the 
elbow- joint axis of rotation. By choice, r 2 = SN = 6 in. in this example. The 
true modified parameters are shown in table II for different values of a 2 (angle 
between rotational axes of the elbow and shoulder joints) to locate the elbow-joint 
axis system at point E in figure 7. Errors in the calculated parameters using 
measurement data rounded off to the nearest ten-thousandth, one-thousandth, and one- 
hundredth of an inch are also shown. In addition to locating the elbow— joint axis 
system at the desired point E in figure 7, the modified parameters can be accu- 
rately calculated even as the lines of rotation of the elbow and shoulder joints 
approach a parallel condition (a 2 0°). 


CONCLUDING REMARKS 

At present, the most popular way to describe the relative location of successive 
joint axis systems in a robot arm is to use the Denavit-Hartenberg parameters. How- 
ever, a recent justifiable criticism is that one of these parameters approaches infi- 
nity when two successive joints have nearly parallel rotational axes. Geometrically, 
this parameter removes the joint axis an excessive distance from the robot arm; com- 
putationally, this large parameter leads to an ill-conditioned transformation matrix. 
In this paper, a simple modification in the location of this axis system easily over- 
comes this criticism. This modification results by insisting that a transverse vec- 
tor between successive joint rotational axes be normal to one of the rotational axes 
instead of to both axes. This simple modification leads to modified Denavit- 
Hartenberg parameters that favorably locate successive joint axis systems. 

An example is given with respect to the rotational axes of the elbow and shoul- 
der joints in a robot arm. The regular and modified Denavit-Hartenberg parameters 
(that locate the elbow- joint axis system relative to the shoulder- joint axis system) 
are extracted by an algebraic method via simulated measurement data. 

A point near the robot hand (off the line of rotation) generates a circular 
trajectory as the elbow joint is rotated. Three position vectors to this point and 
the corresponding three elbow joint angles (referenced to an initial position) are 
simulated for the extraction process. 

For small misalignments of the shoulder and elbow joints away from a parallel 
condition (i.e., parallel rotational axes), the Denavit-Hartenberg parameters locate 
the elbow-joint axis system far away from the robot arm; the modified parameters lo- 
cate the axis system at the desired place on the robot arm. In addition, for a given 
accuracy of the measurements used in the parameter-extraction process, the extracted 
values for the Denavit-Hartenberg parameters yielded considerably larger errors than 
did the extracted values for the modified parameters. It appears that the modified 
parameters provide a more natural location of successive joint axis systems and are 
useful in the industrial calibration of robot arms. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
March 13, 1986 
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APPENDIX 


EQUATIONS FOR SIMULATED MEASUREMENTS 

By Euler 1 s theorem (ref. 5), a vector r rotated by an angle 0 about a line 
(or unit vector ( 5 ) becomes 


r* = cos 0 r + r • w(1 - cos 0)w - (r x w) sin 0 (A1 ) 

This equation is applied to obtain an expression for the location of point W on the 

robot arm in figure A1 with respect to the waist axis system as the elbow joint 9^ 
is rotated for a specified value of cx 2 . 

For the position of the robot arm in figure A1 , and with respect to the X, Y, 

and Z axes, let r = (0, WE, 0) T be a vector from E to W, and let 

5 = = (sin a 2 , 0 , cos a 2 ) T be a vector along the rotational axis of the elbow 

joint. Moreover, let 0 = ©3 be the elbow- joint angle. Then r* gives the new 
location of W as a function of 0g and a 2 . The location of W with respect to 
the waist axis system is then simply 


( A2 ) 


where r^, r^, and are the components of r* (or coordinates of W) along the 

X, Y, and Z axes. 

Simulated measurement vectors A, B, and C are the same as W, except 0 3 is 
replaced by 0 A , © B , and 0 C , respectively. For the example in the text of this 

report, 0 A = 0°, 0 B = 45°, and 0 C = 90°. 
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Figure A1 . - Elbow- joint rotational axis, with 
respect to an auxiliary axis system, for 
simulating measurement data. 
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TABLE I.- ERRORS IN CALCULATED DENAVIT-HARTENBERG PARAMETERS THAT 
LOCATE ORIGIN OF ELBOW-JOINT AXIS SYSTEM FOR ROBOT ARM 

[ Simulated measurements to point on robot arm correspond] 
to elbow- joint angles of 0°, 45 °, and 90° J 


True simulated value of - 

Simulated point 
measurements 
rounded off to 
nearest - 

Absolute error in 
computed value of - 

a 2' deg 

r 2 , in. 

r • 

r 2 / in. 

a2 / m . 

0,01 

-97396.82 

0.00 

10 -4 

11 .22 

0.30 

.01 

-97396.82 

.00 

10 -3 

6410.25 

1 .83 

.01 

-97396.82 

.00 

10 -2 

97406.82 

6.51 

0.10 

-9734.27 

0.00 

10" 4 

9.58 

0.00 

.10 

-9734.27 

.00 

10“ 3 

55.45 

.30 

.10 

-9734.27 

.00 

10 -2 

808.07 

1 .82 

1 .00 

-967.93 

0.00 

10" 4 

0.37 

0.00 

1 .00 

-967.93 

.00 

10 -3 

1 .10 

.00 

1 .00 

-967.93 

.00 

10~ 2 

4.46 

.28 

10.00 

-90.41 

0.00 

10“ 4 

0.00 

0.00 

10.00 

-90.41 

.00 

10~ 3 

.04 

.01 

10.00 

-90.41 

.00 

10~ 2 

.25 

.03 
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TABLE II.- ERRORS IN CALCULATED MODIFIED PARAMETERS THAT LOCATE ORIGIN 
OF ELBOW-JOINT AXIS SYSTEM FOR ROBOT ARM WHEN TRANSVERSE VECTOR 
IS CONSTRAINED NORMAL TO ELBOW-JOINT ROTATIONAL AXIS 

[ Simulated measurements to point on robot arm correspond 1 
to elbow- joint angles of 0°, 45°, and 90° J 


True simulated value of - 

Simulated point 
measurements 
rounded off to 
nearest - 

Absolute error in 
computed value of - 

a 2 , deg 

$2' in - 

r\ 2 , in. 

*29 in. 

l 2 > in * 

^2 1 in • 

0.01 

0.00 

17.00 

6.00 

10* 4 

0.0005 

0.0005 

.01 

.00 

17.00 

6.00 

10 -3 

.0006 

.0006 

.01 

.00 

17.00 

6.00 

10~ 2 

.0028 

.0028 

0.10 

0.00 

17.00 

6.00 

10~ 4 

0.0001 

0.0001 

.10 

.00 

17.00 

6.00 

10” 3 

.0006 

.0006 

.10 

.00 

17.00 

6.00 

10~ 2 

.0028 

.0028 

1 .00 

0.00 

17.00 

6.00 

10~ 4 

0.0002 

0.0002 

1 .00 

.00 

17.00 

6.00 

10~ 3 

.0003 

.0003 

1 .00 

.00 

17.00 

6.00 

10“ 2 

.0041 

.0067 

10.00 

0.00 

17.00 

6.00 

10“ 4 

0.0001 

0.0001 

10.00 

.00 

17.00 

6.00 

10“ 3 

.0004 

.0006 

10.00 

.00 

17.00 

6.00 

10“ 2 

.0022 

.0035 









Figure 2.- Initial position of robot arm and 
joint axis systems. 


23 









(a) 


*> 

D 2 normal to shoulder rotational axis . 




Elbow 

rotation 


0 


2 



Shoulder 

rotation 


\**~ r 2 + 

(h) normal to elbow rotational axis Z 2 . 

Figure 6.- Location of point E from N by h 2# 



H 



Figure 7.- Geometry used in locating 
elbow-joint axis system. 
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